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METHOD FOR CALCULATION OF LAMINAR HEAT TRANSFER IN AIR FLOW AROUND 
CYLINDERS OF ARBITRARY CROSS SECTION (INCLUDING LARGE TEMPERATURE 
DIFFERENCES AND TRANSPIRATION COOLING) : 


By E. R. С. Ескевт and Јонм N. B. LrviNaoop 


SUMMARY 


The solution of heat-transfer problems has become vital for 
many aeronautical applications. The shapes of objects to be 
cooled can often be approximated by cylinders of various cross 
Sections with flow normal to the axis as, for instance, heat 
transfer on gas-turbine blades and on air foils heated for deicing 
purposes.. A laminar region always exists near the stagnation 
point of such objects. 

A method previously presented by E. В. G. Eckert permits 
the calculation of local heat transfer around the periphery of 
cylinders of arbitrary cross section in the laminar region for 
flow of a fluid with constant property values with an accuracy 
sufficient for engineering purposes. The method is based on 
exact solutions of the boundary-layer equations for incompres- 
sible wedge-type flow and on the postulate that at any point on 
the cylinder the boundary-layer growth is the same as that on a 
wedge with comparable flow conditions. This method 18 
extended herein to take into account the influence of large tem- 
perature differences between the cylinder wall and the flow as 
well as the influence of transpiration cooling when the same 
medium as the outside flow is used as coolant. Prepared 
charts make the calculation procedure very rapid. For cylinders 
with solid walls and elliptic cross sections, a comparison 18 
made of the results of calculations based on the presented method, 
the resulis of calculations by other known methods, and МИНЕ 
obtained in experimental investigations. 


INTRO DUCTION 


Calculation of the heat transferred to cylinders with : 


arbitrary cross sections from air flowing normal to the axis 
by a solution of the boundary-layer equations is a difficult 
problem, even when the laminar region is considered. The 
problem is especially complicated by the large number of 
parameters influencing heat transfer. Such parameters are: 
the shape of the cross section of the cylinder, the Mach 
number which determines the flow outside the boundary 
layer, the temperatures on the surface of the cylinder as 
well as in the stream, the stream velocity determining the 
internal heat generation, and the temperature distribution 
around the circumference of the cylinder. If the cylinder 
is cooled by the transpiration-cooling method in which a 
coolant is ejected through a porous surface into the outside 
stream, the amount of coolant and its distribution around 
the circumference of the cross section of the cylinder are 
additional parameters. Even if a solution is obtained for 


, Prandtl number. 


such a problem, for instance by use of an electronic computer, 


the solution is very restricted because of the many parameters. 
Up to the present time, therefore, the problem has been 
attacked only under simplifying restrictions. 

The restrictions most commonly used are: (1) low veloci- 
ties, (2) constant property values, (3) constant wall tem- 
peratures, and (4) impermeable surfaces (no transpiration 
cooling). Under restriction (2), the development of the 
boundary layer along the cylindrical surface is independent 
of the heat transfer; available knowledge on the flow bound- 
ary layer can therefore be used as a basis for a heat-transfer 
calculation. Under the simplifying assumptions, which are 
necessary in order to transform the general viscous-flow 
equations into the boundary-layer equations, the develop- 
ment of the flow boundary layer does not depend immediately 
on the shape of the cross section of the cylinder but only on 
the velocity distribution in the stream outside the boundary 
layer and along its surface. 

One method which was applied successfully to obtain a 
solution of the flow boundary-layer equation developed the 
stream velocity. along the surface of the cylinder in a power 
series of the distance from the stagnation point measured 
along the circumference of the cylinder. In reference 1, this 
method is used to solve the heat-transfer problem. It is also 
shown that the temperature field within the boundary layer 
can be presented in a power series of the distance from the 
stagnation point in which the single terms contain only uni- 
versal functions of a dimensionless wall distance and of the 
Prandtl number of the fluid. The heat transfer to the surface 
is given by an analogous series with terms depending on the 
The calculation of the universal functions, 
however, is а tedious process, and accordingly these functions 
are known only for a limited number of terms. For.air with a 
Prandtl number of 0.7, they are presented inreference 1. For 
& gas with a Prandtl number of 1, they are contained in ref- 
erence 2, which is based on reference 3, in which the boundary- 
layer flow on a yawed cylinder is calculated. The fact 
that the boundary-layer equation for the velocity component 
parallel to the axis of a yawed cylinder is identical in form to 
the boundary-layer equation describing the temperature field 
for a fluid with a Prandtl number of 1, flowing normal to the 
axis of the cylinder, was used in reference 2 to determine 
heat transfer to such cylinders. The presentation of more 
terms of the series is announced in reference 4. It was found, 
however, that the velocity distribution for only a limited 


1 Supersedes МАСА TN 2733, "Method for Calculation of Heat Transfer In Laminar Region of Air Flow Around Cylinders of Arbitrary Oross]Section (Including Large Temperature Differs 


ences and Transpiration Cooling)” by E. R. G. Eckert and John N. B. Livingood, 1952, 
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range of cross sections of cylinders can be represented by а 
power series converging rapidly enough that the number of 
the known universal functions is sufficient to calculate the 
heat transfer. i 

The difficulties connected with a solution of the boundary- 
layer equations point out the need for an approximate 
approach by which, with a small expenditure of time, heat- 
transfer coefficients can be determined with an. accuracy 
sufficient for engineering purposes. A considerable number 
of such approaches have been tried in the past; the results 
differ greatly as shown in figure 1, taken from reference 2. 


Heat-transfer coefficient 


` Flat plote-- 


ША 2 .3 .4 25 6. T 
Distance from stagnation point 
FIGURE 1.—Heat-transfer coefficients for cylinder (ref. 2). 


The simplest procedure is probably that in which the heat- 
transfer coefficients as calculated in reference 5 are used for 
flow of constant velocity along a flat plate. The fact that in 
reality the stream velocity varies along the cross section of the 
cylinder is taken into account by calculating the local heat- 
transfer coefficients by use of the velocity found in the stream 
at the considered distance from the stagnation point. This 
method is contained in a summary presented in reference 6. 
Unfortunately, such an approach gives heat-transfer coeffi- 
cients which are considerably low in many cases (see fig. 1). 

Better agreement was obtained by another approach (ref. 
7) which uses, instead of the flat-plate solution, & family of 
solutions of the boundary-layer equations which can be 
obtained in a general form, namely, for the case where the 
stream velocity varies along the surface as a certain power of 
the distance from the stagnation point. Such a velocity 
variation is obtained in incompressible flow around wedges. 
The solutions for such a type of flow were used to obtain 
approximate heat-transfer coefficients for a cylinder with 


arbitrary cross section by stipulating that the local heat- 


REPORT 1118—NATIONAL ADVISORY COMMITTEE FOR AERONAUTICS 


transfer coefficient on any location along the cylinder is 
identical with the local heat-transfer coefficient on & wedge 
for which, at the same distance from the stagnation point, 
the stream velocity and its gradient are the same as those on 
the investigated cylinder. This approach was subsequently 
used by different authors, and is described, for instance, in 
references 8 and 9. It takes into account the:stream condi- 
tions which influence the boundary-layer growth at the 
location at which the heat transfer is going to be determined; 
however, it does not properly account for the development of 
the boundary layer in the range upstream of the point con- 
sidered. This development may be different on the cylinder 
and on the equivalent wedge. 

Another group uses an integrated momentum equation 
for the boundary-layer flow as proposed by von Kármán and 
K. Pohlhausen (refs. 10 and 11, respectively) to calculato 
the velocity boundary layer. Different procedures were 
proposed for determining local heat-transfer coefficients 
from the known velocity boundary layer. Some investi- 
gators use Reynolds analogy directly (ref. 12) or with a 
correction for Prandtl numbers different from 1 (ref. 13). 
Such approaches give heat-transfer coefficients which aro 
considerably high in many cases, as shown, for instance, in 
figure 1. More accurate results were obtained when the 
heat transfer was determined by solving an integrated heat- 
flow equation for the boundary layer. The velocity field 
within the boundary layer has to be known in this approach, 
since the flow velocities within the boundary layer occur in 
the mentioned heat-flow equation. This method was origi- 
nated by Kroujiline (ref. 14). Extensions and simplifica- 
tions are contained in references 15 to 18, and an extension 
to compressible flow of & fluid having & Prandtl number 
equal to 1 is found in references 19 and 20. Useful informa- 
tion is also contained in a summarizing report (ref. 21). 

Another approach is based on the fact that the use of tho 
heat-transfer coefficients for wedge-type profiles as described 
previously was found to give fairly accurate heat-transfer 
coefficients. “It should be expected that these heat-transfer 
coefficients can be improved to a degree whith is sufficient 
for all engineering purposes by a method which takes into 
account in some approximate way the previous history of 
the boundary layer. Such a method, called the equivalent 
wedge-type flow method, is proposed in reference 22, ex- 
tended to heat transfer at high flow velocities and variable 
wall temperature in reference 23, and extended to transpira- 
tion cooling with small temperature differences in reference 
24. The advantages of this method are that no knowledge 
of the velocity boundary layer is required and that it can 
be readily extended to take into account the effects of largo 
temperature differences, of transpiration cooling, and of 
variable wall temperature as soon as the corresponding 
solutions for the wedge-type flow are available. 

Such an extension was made at the NACA Lewis labora- 
tory during 1950-51 and is described herein. It is based on 
exact boundary-layer solutions for wedge-type flow with 
large temperature differences and with transpiration cooling 
(refs. 25 and 26). Charts were prepared which make the 
calculation of heat transfer around cylinders of any arbitrary 
cross section more rapid. 
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SOLUTION OF BOUNDARY-LAYER EQUATIONS FOR WEDGE- 
TYPE FLOW 


BOUNDARY-LAYER EQUATIONS 


The following boundary-layer equations describe the 
velocity and temperature fields in a laminar steady two- 
dimensional gas flow: the momentum equation, the conti- 
nuity equation, and the energy equation. The momentum 
equation is 


ди д oe _ Op 
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when body forces are neglected. (All symbols are defined 
in appendix A; consistent units are used throughout the 
report. Since the pressure variation normal to the surface 
throughout the boundary layer may be neglected, it follows 
that the pressure is prescribed by the conditions in the stream 
outside the boundary layer and can be connected with the 
velocity u, in the stream and just outside the boundary 
layer by the Bernoulli equation 
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The introduction of this expression changes Ше momentum 
equation to the form 
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The continuity equation is 
д д . | 
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and the energy equation is 
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The heat. generated by internal friction, described by the 
second term on the right side of equation (4), and the temper- 
ature variation connected with expansion, described by the 
third term, can be neglected as long as the difference between 
the total and the static temperature in the gas stream is small 
compared with the difference between the wall temperature 
and the temperature in the gas stream. For this condition, 
then, only the first term on the right side of equation (4) is 
retained, and the energy equation assumes the form 


pe; (и Sete 92) 5. (655 (5) 


Equations (2), (3), and (5) include the case of transpiration 
cooling when the same medium as that in the outside flow is 


used as coolant and the boundary conditions are properly 
defined. 
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The property values и, Е, ср, and р appearing in the equation 
depend on temperature and pressure. The variation with 
pressure can be neglected at the low velocities to which the 
energy equation was already restricted by disregarding the 
internal friction and the expansion terms. The influence of 
the temperature dependency, however, may be appreciable 
in applications with large temperature differences within the 
boundary làyer. Solutions of the boundary-layer equations 
which take into account the temperature variation of the 
property values were obtained in references 9, 25, and 26, in 
which the partial differential equations were transformed into 
total differential equations. 


CHANGE OF VARIABLES 


The transformation of the partial differential equations 
into total differential equations is possible under the following 
specialized conditions: The stream velocity is assumed to 
vary as a power function of the distance from the stagnation 
point measured along the surface of the cylinder. 


It has recently become customary to refer to the exponent m 
in this equation as “Euler number." The Euler number сап 
be expressed by the Bernoulli equation in the following way: 


_әр= . 


pus] x 


(8) 


In addition, the temperature of the wall is assumed to be 
constant and the property values are assumed to vary pro- 
portionally to a power of the absolute temperature T. The 
numerical calculations were made for air. 'The exponents 
used were 0.7 for the viscosity, 0.85 for the heat conductivity, 
0.19 for the specific heat, and — 1.0 for the density. 

The variables 
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are used to transform equations (2), (3), and (5) into total 
differential equations presenting f and 0 as functions of у only. 
The stream function y appearing in equations (9) is defined in 
such а way as to eliminate the continuity equation (3). 


u= Aen 
(10) 
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Introducing the new variables into the second of equations 
(10) leads to the following expression for the velocity com- 
ponent normal to the surface: 


тезі? Pur fm—10*?| (11) 


--- +e Se eee eer 


226 


The velocity at the surface itself follows: 


Itm pr ae (12) 
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The transformation therefore e & certain variation 
of the coolant velocity v, along the surface, since the function. 
fo has to be constant (independent of х). The stream 
velocity is described by equation (7); thus, the coolant 
velocity vo, is also proportional to some power of х. Such a 
variation of the coolant velocity leads to a constant wall 
temperature and is therefore consistent with the assumed 
constant wall temperature when heat transfer by radiation 
may be neglected (ref. 27). The transformed equations are 
presented in references 9, 25, and 26, together with the solu- 
tions for a Prandtl number Pr of 0.7, and for a range of 
Euler number m, temperature ratio Т,/Т,, and the parameter 
Ж describing the cooling-air flow through a porous surface. 
The results contain expressions for the thickness of the flow 
boundary layer which are defined in two ways: the displace- 


ment thickness 
im L (1-2 u (13) 
and the momentum thickness 
= 0 rare 841 = dy (14) 


The thermal boundary layer is characterized in this report 
by the convection thickness 


ри =, Ба dy 
` Jo реш, Т,-1 — T; 


In addition, a thermal boundary-layer thickness, which is 
defined as follows, will be used herein: 


=], (naa) 


Values for this boundary-layer thickness can be easily calcu- 
lated from results presented in references 25 and 20. 
APPLICATION TO HIGH VELOCITIES 
The solutions described apply exactly only to flow with 
low velocities. Practically, the limiting velocity up to which 
it is possible to neglect the frictional and the expansion terms 
can be set quite high for a gas; this fact can be understood 


(16) 


from the following transformation of the energy equation,. 


in which only the specific heat is regarded constant. If the 
momentum equation (1) is multiplied by the velocity и and 
added to the energy equation (4) and if, in addition, the total 
temperature T'p=7'+-u?/2c, is introduced, the following ex- 


pression is obtained: 
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The last term on the right side of the equation vanishes for а 
Prandtl number equal to 1. In this case, the energy equa- 
tion has the same form as the one for low velocities in which 
the friction and the expansion terms were neglected. The 
only difference lies in the fact that the totol temperature 
appears in the energy equation. When the Prandtl number 
is approximately 1, the last term in equation (17) will be 
comparatively smell up to considerable velocities, and the 
energy equation (5) used in the following considerations 
applies to this condition when the temperature T is inter- 
preted as total temperature. It will be shown later that as 
far as heat transfer is concerned, the range in which the 
results of & calculation with equation (5) may be used can 
be extended even further by using & properly defined adia- 
batic wall temperature instead of the total gas temperature. 


The property values и, k, ср, and p depend for gases on the 
temperature. 'This dependency was taken into account in 
the -described calculations. The density depends, in addi- 
tion, on the pressure, and the pressure variation may become 
considerable at high Mach numbers. There are indications, 
however, that calculations which neglect this pressure varia- 
tion can bé used with sufficient accuracy over the entire 
subsonic range, as is pointed out in reference 28, in which 
an investigation of results obtained by L. Howarth (ref. 29) 
is reported. 


EXTENSION OF THEORY TO ARBITRARY BODIES 
DETERMINATION OF EQUIVALENT WEDGE 


The solutions discussed in the previous paragraph are in an 
exact sense restricted to a-certain type of velocity variation 
along a cylindrical surface, namely, a stream velocity which 
just outside the boundary layer is proportional to some power 
of the distance from the stagnation point. Such a velocity 
distribution is realized, for instance, in incompressible flow 
around wedges. The wedge-type solutions may be used, 
however, to obtain approximate heat-transfer coefficients on 
cylinders of arbitrary cross section. In one approach in this 
direction, it is assumed that the heat-transfer coefficient on 
any point along the circumference of & profile with arbitrary 
cross section is the same as that on a wedge at the same dis- 
tance from the leading edge, provided the stream velocity and 
its gradient on the wedge and on the arbitrary profile have 
the same value at the location considered and that the tem- 
perature ratio T,/T,is Ше same. It will be shown that such 
an approach takes into account the right stream conditions 
at the local spot for which the heat-transfer coefficient is to 
be determined. However, the previous history within the 
boundary layer is not properly considered. Numerical 
calculations presented herein show that heat-transfer co- 
efficients obtained in such a manner are in most cases within 
about 15-percent agreement with experimental data. It is 
to be expected that a modification which accounts in some 
approximate manner for the conditions in the boundary layer 
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upstream of the point under consideration should improve 
this approximation to the desired degree. This modification 
is made in reference 22 by the stipulation that the rate of in- 
crease of the boundary-layer thickness is the same on the 
considered point of an arbitrary profile and on the point of a 
wedge which has the same boundary-layer thickness, the 
same stream velocity, and the same stream velocity gradient. 


This same stipulation will be used in the present report. For- 


a given temperature ratio Т./Т., the heat-transfer coefficients 
on a wedge depend ор the Euler number m and the value ў, 
characterizing the coolant flow through & porous surface. 
These parameters which define the equivalent wedge profile 
will now be expressed by the boundary-layer thickness and 
the local stream velocity gradient. 


For the wedge-type profile, the stream velocity is expressed 
by the power law 
=O" (18) 


in which the value £ expresses the distance from the leading 
edge measured along the wedge surface in order to distinguish 
it from the distance of the point under consideration from 
the stagnation point on the arbitrary profile, which is denoted 
by z. The variables used for the transformation of the 
original boundary-layer equations in the previous section 
may now be written 


€ DU 
n=y y p (19) 
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Corresponding to а certain value y, which indicates the 
boundary-layer thickness 6, there is a value mof the coordinate 
7 defined by the equation 


пь--6 " Pwlls 
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In order to eliminate the distance £ from this equation, equa- 
tion (18) is differentiated to obtain 


(21) 
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Since the velocity gradient on the wedge profile is assumed 
the same as that on the profile under consideration, it follows 
that ди, = 0и, |02. This eee gives for the coordinate 
£ the expression 
оти 
Е даа (43) 


When this expression is introduced into equation. (21), there 


is obtained 
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In this expression, m (denoted as (5/z)/Re in refs. 25 and 26) 
is & function of the Euler number m and of the coolant-flow 


parameter Ju. Therefore, if this equation is written in the 
form 
2, Pd dus ' (9 
Път "m dz (24) 


the left side is a function of m and fo, and equation (24) re- 
lates both values to the boundary-layer thickness 8 and the 
velocity gradient du,/dz. In order to obtain a second relation 
for m and fw, the coordinate Е is replaced in equation (20). 
The result is 


oni a (25) 
which is written again in such a way that the left side is a 


function of the Euler number m and the flow parameter f, 
which can be calculated from the results in references 25 


and 26. Both equations (24) and (25) are therefore suffi- 


cient to determine the equivalent wedge profile. 


EQUATIONS.FOR BOUNDARY-LAYER THICKNESS AND HEAT TRANSFER 


The next step is to develop a differential equation for the 
boundary-layer thickness from the postulate that the 
boundary-layer gradient 48/45 is the same for the real profile 
as for the equivalent wedge profile. For the wedge profile, 
the boundary-layer thickness is given by the expression 


M 
— Це B 
= Орь 


which is obtained from equation (21) by replacing the stream 
velocity with equation (18) and solving for the boundary- 
layer thickness. A differentiation of this equation and the 
use of equations (23) and (24) result in 


dz 2m Us а (26) 


° р.и, 
This is a differential equation for the boundary-layer thick- 
ness which contains only values which are known for the 
profile under consideration or which are determined from 
equations (24) and (25) for the equivalent wedge-type flow. 
An integration of the differential equation gives the boundary 
layer along the circumference of the profile under consider- 
ation. 

The local heat-transfer coefficient is defined by the follow- 


ing equation: 


Introducing the dimensionless temperature ratio given in 
equation (9) and the coordinate £ results in 
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The heat-transfer coefficient may be calculated from this 
expression as soon as Ше boundary-layer thickness ô is 
known, since 6, and m are functions of m and f, contained in 
references 25 and 26. 

Up to the present time no recommendation has been 
made as to which boundary-layer thickness should be used 
in the prescribed procedure. When the momentum thickness 
is used in the foregoing equations, it is easily understandable 
that the integrated momentum equation is satisfied and 
that the method of calculation becomes the same as the one 
proposed Бу von Kármán in reference 10. This fact can be 
proved mathematically by a procedure completely analogous 
to the one used in appendix B. On the other hand, the use 
of the convection thickness as defined in equation (15) 
satisfies the integrated heat-flow equation within the bound- 
ary layer, аз shown in appendix B. The use of both boundary- 
layer thicknesses leads to somewhat different results for the 
local heat-transfer coefficient, and a question arises as to 
which is preferable. It is pointed out by Schuh in reference 
23 that, for the purpose of determinmg heat-transfer co- 
efficients, it is more important to satisfy the heat-flow 
balance; the use of the convection thickness was therefore 
recommended. In reference 22, the use of the thermal 
boundary-layer thickness as defined in equation (16) is 
investigated, and the results of the calculation with this 
boundary-layer thickness are found to agree even better 
with measured values and with other calculations. The 
convection thickness 6, and the thermal thickness 8, for the 
boundary layer will therefore be used in pone in the 
following numerical evaluations. 


CALCULATION PROCEDURE 
USE OF DIMENSIONLESS VARIABLES 


The procedure which may be followed in determining 
local heat-transfer coefficients with the relations developed 
in the preceding section is now explained. Figure 2 shows a 


7s 


Stagnation / ` 
point - 


FIGURE 2.—Sketch of cylinder indicating notation used, 


sketch of a cylinder with arbitrary cross section and the 
notation used in the analysis. Before numerical calculations 
are made, however, it is advisable to change to dimensionless 
quantities. ln order to make this change, the distance z is 
divided by the major axis L of the cylinder and the mass 
velocity in the direction of x is divided by an upstream mass 
velocity. All lengths and mass velocities parallel to y are, 
in addition, multiplied by the square root of the Reynolds 
number Re, based on the major axis and the upstream mass 
velocity. The dimensionless variables which are sub- 
sequently needed are 
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By use of these dimensionless quantities, equation (26) is 
transformed into 
а“ M 


da* Sut 125) 
where 
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according to equations (24) and (25), which, in dimensionless 
values, are 
ди, ще 
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and 
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Introduction of the dimensionless quantities into equation 
(27) leads to 
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and 
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CHARTS AND CALCULATION PROCEDURE FOR PRESCRIBED COOLANT 
FLOW 

Charts have been prepared which present the functions 
M and N as expressed by equations (34) and (37) in depend- 
ence on (Фи /дт“)6“ and 289“. The charts presented herein 
were constructed from results presented in references 25 
and 26. In figures 3 and 4, the dimensionless convection 
thickness of the boundary layer is used; in figures 5 and 6, 
the dimensionless thermal boundary-layer thickness is used. 

At the stagnation point of any blunt nosed cylindrical 
body, conditions are the same as those at the stagnation 
point of a plate normal to the flow. Therefore m=1, but 
the value of 6* is unknown. However, there exists at the 
stagnation point a unique relation v$5*— F((du?[dz*)5*?] 
which may, for instance, be read along the abscissa in figure 
3 or in figure 5. Squaring this equation and dividing both 
sides by (du? /dz*)5* result in 


duf|dz* (ди ааа Када" 
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dus «о 
(Ы) ТУТ», 2. dx 8 
FIGURE 3.—Chart for use in determination of M for dimensionless convection boundary-layer thickness. Pr, 0.7. 
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(с) TH Tu, 4. 
FIGURE 3.—Concluded. Chart for use in determination of M for dimensionless convection boundary-layer thickness. Pr, 0.7. 
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(a) E 55 1. 
FIGURE 4.—Chart for use in determination of N for dimensionless convection boundary-layer thickness. Pr, 0.7. 
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FIGURE 4.—Concluded, Chart for use in determination of N for dimensionless convection boundary-layer thickness. Pr, 0.7. 


(а) TH Te, 1. 
FIGURE 5.—Chart fof use in determination of M for dimensionless thermal boundary-layer thickness. Pr, 0.7. 


These relations are presented in figure 7 for the dimension- 
less convection boundary-layer thickness and in figure 8 for 
the dimensionless thermal boundary-layer thickness. 

By use of these charts, the calculations for any profile 
can be made in a very simple manner for either the dimen- 
sionless convection or the dimensionless thermal boundary- 
layer thickness. The method of solution for the convection 
thicknesses is described subsequently. For the thermal 
thickness, the procedure is the same. 

The values-of u, and du,/dz must be found for the cylinder 


profile under consideration either by measurement or by a 
solution of the inviscid-flow equations. The coolant veloc- 
ity 0, is prescribed by the porosity of the wall and by the 
pressure distribution around the profile. From these terms, 
the values of uf, duf/dx*, and vi can be calculated. The 
value of ôf at the stagnation point can be determined from 
figure 7 in the following way: The value of vz^/(du7 ах“) is 
computed, and the corresponding value of (duf/dz*)8?^ is 
read from figure 7. A simple algebraic operation then yields 
the desired value of 6 at the stagnation point. 
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FIGURE 5.—Continued. Chart for use in determination of M for dimensionless thermal boundary-layer thickness. Pr, 0.7. ` 


The dimensionless convection boundary-layer thickness 
б; along the cylindrical surface is determined from equation 
(33); for the numerical evaluation presented herein, this 
equation was solved by the method of isoclines with the aid 
of figure 3, depending upon which ratio of stream to wall 
temperature is applied. Equation (33) determines the direc- 
tion of the tangents to the different 57-curves which satisfy 
the equation. The task is to find that curve which contains 
the &-value previously calculated for the stagnation point. 
Гог chosen values of z* and 62, values of (du?/dx*)8? and 
zê, are computed and the value of M is read from the 
appropriate part of figure 3. Equation (33) then gives the 
slope of the tangent at this selected value of 2“ for the 
assumed ôf. Several values of 6* are used for this z*. The 
same calculations are repeated for other values of z*. If 
the chosen distance between these z*-values is small enough, 
an accurate curve of à; against z* can be drawn which starts 
at the desired previously calculated value of 57 at the stag- 


nation point and which will have the correct slope at each 
value of z* considered. Figure 9 illustrates this method of 
solution. "Values of N can then be obtained for each of the 
correct à;-values and the considered v*-value for each z* 
from figure 4 after (du7/dz*)5;? and vi; are computed (the 
ratio of stream to wall temperature under consideration 
determines which part of figure 4 should be used). "The 


value of Nu Re, can finally be obtained from equation (35). 


The same calculation procedure can be used when the 
dimensionless thermal beundary-layer thickness is considered. 
Figure 8 is used for the determination of the value of 5* at the 
stagnation point; figure 5 is used to determine M; and figure 
6 is used to determine N. The particular ratio of stream to 
wall temperature under consideration determines which parts 
of these figures apply for the calculation of the values of AZ 
and М. Finally, equation (35) gives the desired value of 
Nu Reo. 
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FIGURE 5.—Concluded. Chart for use in determination of M for dimensionless thermal boundary-layer thickness. Pr, 0.7. 
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FIGURE 6.—Chart for use in determination of N for dimensionless thermal boundary-layer thickness. Pr, 0.7. 
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FIGURE 6.—Concluded. Chart for use In determination of N for dimensionless thermal boundary-layer thickness, Рт, 0.7. 
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FIGURE 7.—Chart for use in determination of dimensionless convection boundary-layer thickness at stagnation point. 
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FIGURE 8.—Ohart for use in determination of dimensionless thermal boundary-layer thickness at stagnation point. 
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FIGURE 9.—Isocline solution of boundary-layer equation. 


CHARTS AND CALCULATION PROCEDURE FOR PRESCRIBED WALL 
TEMPERATURE 


'The heat-transfer coefficients determined by the values of 


Nu Rey can now be used to calculate the surface tempera- 
ture of the cylinder when the outside stream temperature 
and the temperature at which the coolant is supplied to 
the interior of the cylinder are known. For this purpose, & 
heat balance for an element of the wall as shown in figure 10 
is set up. The cylindrical volume element considered may 
have two plane surfaces, one surface (1) coinciding with the 
‘outside surface of the cylinder wall and the other (2) apart 
from the inside surface of the wall by such a distance that it 
is situated outside the boundary layer present on this side. 
(The inside surface has to be considered as a surface of a 
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FIGURE 10.—Cross section through part of cylinder wall used in setting up heat balance. 
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wall to which suction is applied and on which а, boundary 
layer builds up as shown in ref. 30.) The mantle surface 
(3) of the cylinder may be normal to the wall surfaces. 
Heat is carried by convection with the cooling air through 
surfaces 1 and 2. The amount per unit time is indicated in 
figure 10. It is assumed that the coolant is heated up to 
the wall surface temperature Т, when it leaves the wall. 
This assumption is usually well fulfilled. Heat will also be 
transferred by conduction through the fluid layers immedi- 
ately adjacent to the outside wall surfaces, the amount being 
—k,, (ӘТ/оу).а А. In addition, heat may be transferred to 
the outside wall by radiation; it may be q, dA. Heat may 
also flow into the volume element by conduction in the solid 
material or by transverse flow of the cooling air., The sum 
of all these individual flows may be 4: dA. Then the heat 
balance is 


| DIN — 
q:d- Pet Cp Po Tc— Cp Duo T w— Ка ЕЗ 


The heat —k, (07'/Oy). transferred per unit area from the 
gas to the wall is expressed in this report by а heat-transfer 


coefficient 
> 
МТ,-Туес-і, (= ) 


Combining these two equations results in 


91-9 A(T,— To) 65 petu(Tw— Т.) (38) 


This equation permits a calculation of the wall temperature 
for any place on the cylindrical surface when the coolant 
velocity v, is prescribed, when the local radiative heat flow 
q, and the conductive heat flow g, are known, and when the 
heat-transfer coefficient № has been obtained. The conduc- 
tive heat flow qe is usually small and can be neglected. Such 
a calculation results in a wall surface temperature which 
generally will vary along the circumference of the cylinder. 
When the variations are large, the temperature distribution 
obtained can be regarded only as an approximation, since 
the wedge solutions (refs. 22, 25, and 26) on which the 
method in this paper is based were obtained for the case of a 


‘constant wall temperature. 


Usually, however, the problem which faces the designer in 
an application is somewhat different from the one treated. 
The purpose of transpiration cooling is mostly to keep the 
wall temperature of some structural element below the limits 
which the material can withstand. On the other hand, the 
amount of coolant almost always must be kept small, which 
means that local overcooling should be avoided. For the 
wall surface under these conditions, a temperature is pre- 
scribed which should be uniform about the circumference of 
the cylinder and the problem is to find that distribution of the 
coolant velocity v, which results in the desired wall tempora- 
ture. Generally, such an investigation requires a trial-and- 
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error procedure which is very involved. The procedure 
becomes simple and straightforward, however, when the 
radiative heat flow q, and the conductive heat flow qe can be 
neglected. Such a solution, then, is also useful as a starting 
point for the trial-and-error procedure when radiation is 
present. 

The heat balance (eq. (38)) can be transformed to 


h =: роби Т.- T, 
Cp pss, 0 Dass, 0 Т,--Т, 


(99) 


when q,=q.=0. The ratio of temperature differences in this 
equation is now a prescribed value. А similar ratio 
(Т,— Т.)/(Т.- T.) often appears in turbine-cooling work and 
is denoted by e. Introduction of this value and conversion 
to dimensionless values result in 


(40) 


N = 
ат»: 


Another expression for Nuj Re ёо is given by equation (35). 
Combining both equations gives 


N= v0%*5* Pr 1—e 


This equation expresses a relation between the parameters 
N and 526“ in figures 4 and 6 which may be used to insert 
lines of constant e into these figures. With the use of these 
lines, the calculation procedure for any specific problem 
becomes quite simple. The procedure will be described for 
Т/То-< 1 (or near 1) and with the use of the convection 
boundary-layer thickness ô. The prescribed temperatures 
fix the value of g. 

At the stagnation point, m=1 and du?/dx* is known. In 
figure 4 (a) the intersection between the line m—1 and the 
line for the prescribed о determines тб, and (duj[dr*)5", 
and, from both values, ô and v? may be calculated. 

The method of isoclines may again be used to determine the 
development of the boundary layer along the cylindrical 
surface. The use of this method implies that the gradient 
dà;[dz* has to be determined for any pair of values 2% and 
ôf. For an assumed 67, the value 99% can be found in 
figure 4 (а) as the value on the prescribed ф-сигуе above the 
known abscissa value (ачғ/ах%) 633, 
М and equation (33), the gradient dó?/dx*. А plot, similar 
to figure 9 determines the boundary-layer thickness, and the 
values vs belonging to these boundary-layer thicknesses 
represent the coolant-flow distribution for the particular 
temperature-difference ratio ¢. 


NUMERICAL EVALUATIONS AND COMPARISONS WITH 
KNOWN RESULTS 
SOLID SURFACES 


The results of the outlined procedure for calculating local 
heat-transfer coefficients have to be compared with experi- 


(41). 


Figure 3 (a) then gives. 
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mental results or calculations by some other method in 
order to check the accuracy. The only cylindrical shape for 
which experimental data or solutions of the boundary-layer 
equations suitable for such a comparison are available seems 
to be the cylinder with & circular cross section. Accordingly, 
local heat-transfer coefficients were calculated according to 
the method proposed in this report by use of the dimension- 
less thermal boundary-layer thickness as well as of the dimen- 
sionless convection boundary-layer thickness. The results 
of these calculations are plotted in figure 11 over the dimen- 
sionless distance from the stagnation point. -Also inserted 
in the figure is a curve representing the average curve through 
the experimentally determined local heat-transfer coefficients 
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FIGURE 11.—Comparison of present method with jar used methods for calculation 
of local heat-transfer coefficients around circular cylinder with impermeable wall p;,0 
Pr, 0. 7; Tu T v, 1. 


mentioned in reference 31. It is shown in reference 22 that 
the measurements correlated well into & single curve when 
the experiments with Reynolds numbers near the critical 
value for transition to turbulence within the boundary layer 
were excluded. The tests with high Reynolds numbers gave 
values of Nujf ео which over the whole upstream side of 
the cylinder were about 10 percent higher then the ones for 
the lower Reynolds numbers. The same behavior is reported 
in гегегепсев 32 and 33 in which it is shown that an increase 
up to 50 percent in the heat-transfer coefficients over the 
expected laminar values was caused by the turbulence level 
in the wind tunnels used. The result of a solution of the 
boundary-layer equation as presented in reference 1 is also 
included in figure 11. 
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This method solves the boundary-layer equations and 
obtains results as a series in the distance along the surface. 
Also inserted are values obtained by use of the Pohlhausen 
flat-plate solution when the free-stream velocity is based on 
the local values and results obtained by the methods of 
references 12 and 13. Heat-transfer coefficients on wedges 
with the same local stream velocity and velocity gradient 
at the same distance from the stagnation point are also 
included (ref. 8). Appendix C explains how these wedge 
solutions were obtained. | 

On а cylinder with a circular cross section, separation 
occurs in the subcritical range near the value z*—0.7. The 
stream velocity distribution around the surface of the cylin- 
der which was needed for the calculations was obtained from 
pressure distributions given in reference 31 and is contained 
in reference 22. 

It may be seen from figure 11 that the use of flat-plate 
values results in heat-transfer coefficients which are consid- 
erably lower than experimental values, whereas the methods 
in references 12 and 13 result in values which are too high. 
Much better agreement is found between the wedge heat- 
transfer coefficients and the experimental results, especially 
near the stagnation point. Farther downstream, the accu- 
racy is improved by the method of this report. For the 
largest distance from the stagnation point, the use of the 
dimensionless thermal boundary-layer thickness results in 
yalues which are higher and the use of the dimensionless 
convection thickness, in values which are lower than the 
experimental ones. The values calculated by Fréssling’s 
solution of the boundary-layer equations are also higher 
than the experimental ones. Fréssling’s method has to be 
considered as an exact solution of the boundary-layer equa- 
tions. In reference 22 it is recommended, on the basis of 
the good agreement between Fróssling's curve and the values 
obtained by the use of the thermal boundary-layer thick- 
ness, that the method of the equivalent wedge flow be based 
on the thermal boundary-layer thickness. The values of 
the heat-transfer coefficients depend primarily on the veloc- 
ity distribution in the stream around the cross section of 
the cylinder. The velocity distribution used for the calcu- 
lation on the circular cylinder is also shown on figure 11. 
The calculations are made for a Prandtl number of 0.7, for 
a solid surface (0,—0) and a temperature ratio T,/T of 1, 
equivalent to the assumption of constant property values. 
These calculations agree within 5 percent with the exact 
calculation and within 8 percent with experiment when the 
immediate neighborhood of the separation point is excluded. 
Similar comparisons have already been made in reference 2 
for a gas with a Prandtl number of 1 and a different velocity 
distribution (see fig. 1). This comparison shows that the 
method proposed by Squire (ref. 16) gives heat-transfer 
coefficients which agree with the exact boundary-layer solu- 
tion to about the same degree as those of the method of 
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the equivalent wedge flow. The same fact holds for the 


‘method indicated in references 15 and 17 especially with 


the improvement given in reference 4. It can be stated m 
summary, therefore, that & number of methods exist today 
which, at least for the circular cylinder, permit the deter- 
mination of heat-transfer coefficients on solid surfaces in the 
laminar region of а gas having constant property values 
with a very good accuracy. The advantage of the equiva- 
lent wedge flow method over those methods just discussed 
is that it gives solutions in а very short time and that it 
can be readily extended to include variable property values 
and transpiration cooling, as was done in this report. The 
wedge solution, according to reference 7, is still more rapid; 
however, the results differ from the experimental values,up 
to 15 percent. 

. Figure 12 gives the analogous results for an elliptic cylinder 
with the axis ratio 1:2. It may be observed that heat- 
transfer coefficients on wedges differ only slightly from those 
obtained for equivalent wedge-type flow, whereas the flat- 
plate values and the ones calculated with references 12 and 
13 are considerably different. No experimental results or 
solutions of the boundary-layer equations for a cylinder with 
such & cross section which could be compared with the ap- 
proximate solutions are known to the authors. The calcu- 
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FIGURE 12.—Comparison of present method with’ formerly used methods for calculation 
of local heat-transfer coefficients around elliptic cylinder with axis ratio of 1:2. v, 0; 
Pr, 0.7; Ta Tu, 1. 
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FIGURE 13,—Comparison of present method with formerly used methods for calculation 
of local DE MN coefficients around elliptic cylinder with axis ratio of 1:4. 05, 0; 
Pr, 0.7; TITL. 


lation with Kroujiline's method, presented in reference 22, 
agrees well with the solutions obtained with the equivalent 
wedge-type flow method. Separation of the flow occurs on 
such a profile near z*—0.8. The stream velocities used are 
calculated values contained in reference 22. 

The agreement between the wedge solutions and the results 
obtained by the method herein is still closer for the elliptic 
cylinder with axis ratio 1:4 (fig. 13). The reason for this 
fact is the type of stream velocity variation occurring on 
elliptic cylinders. Flow separation occurs on this cylinder 
near з#--0.85. The curves in figures 12 and 13 show that 
the stream velocity is comparatively constant over a con- 
siderable part of its circumference after a steep increase near 
the stagnation point. This behavior is more pronounced for 
an axis ratio of 1:4 than for one of 1:2. Ап inspection of 
figure 13 shows that, apart from the region near the stagna- 
tion point, even the flat-plate values give a reasonably good 
approximation. Calculations obtained by use of the dimen- 
sionless thermal boundary-layer thickness extended to the 
flow separation point, whereas those for’ the dimensionless 
convection boundary-layer thickness did not. It therefore 
appears advisable to use the dimensionless thermal boundary- 
layer thickness. | 

Experimental heat-transfer coefficients found at the Uni- 
versity of California for an elliptic cylinder with an axis 
ratio of 1:4 (ref. 34) are about 50 percent higher than the 
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theoretical values shown in figure 13. There are several 
reasons for this discrepancy.. The measured stream-velocity 
distribution was different from the one on which the present 
calculations are based, probably because of a limited width 
of the wind tunnel. The cylinder in the experimental in- 
vestigation was heated by an electric resistance which pro- 
duced a constant heat flow through the surface per unit area. 
Accordingly, the surface temperature varied along the cir- 
cumference of the cylinder, being lowest at the forward 
stagnation point and increasing in the downstream direction. 
Calculations in reference 34 indicate that the higher values 
found in the tests are mostly due to this fact. Another 
increase of the experimental heat-transfer coefficients may 
again be connected with the turbulence level in the wind 
tunnel used, as discussed in connection with the experimental 
results for circular cylinders. 

From figures 11 to 13, it may be concluded that, for 
cylinders with a stream velocity which is fairly constant over 
the greater part of the circumference, local heat-transfer 
coefficients may be obtained with good accuracy from wedge 
solutions. In the region in which the stream velocity varia- 
tion is considerable, the method of the equivalent wedge flow 
gives heat-transfer coefficients with an accuracy sufficient for 


engineering purposes. 
POROUS SURFACES 


Heat-transfer coefficients were calculated by the method 
of the equivalent wedge flow for cylinders with circular and 
elliptic cross sections for transpiration-cooled surfaces and 
different temperature ratios 7,/ 7, by using either the thermal 
or the convection boundary-layer thickness (figs. 14 to 18). 
In these figures it was more expedient to base the Reynolds 
numbers appearing on the ordinate and in the coolant flow 
parameter on the density at wall temperature (Reg) rather 
than on the upstream density (Reo). In figures 11 to 13 and 
19, both Reynolds numbers are identical since they are 
calculated for а temperature ratio 7,/T,—1. The use of 
both boundary-layer thicknesses gives different results only 
for large distances. from the stagnation point. The. vari- 
ation of the heat-transfer coefficients with the ratio of stream 
to wall temperature is comparatively small for solid surfaces. 
This result is in agreement with previous findings. For 
transpiration-cooled surfaces, however, the effect of the 
temperature ratio on the heat-transfer coefficients becomes 
more pronounced, especially on cylinders with nearly circular 
cross sections. In reference 24, the case of transpiration 
cooling with small temperature differences is calculated; this 
reference includes the effect of the temperature ratio by a 
correction factor which is based on the assumption that this 
effect is the same as that determined experimentally for 
impermeable surfaces. A comparison of results shows that 
the procedure in reference 24 underestimates the effect of 
temperature ratio for transpiration-cooled surfaces. In 
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addition, it can be observed that transpiration cooling 


results in a considerable decrease of the heat-transfer co- 


efficients. A larger amount of coolant flow is necessary to 
reduce the heat-transfer coefficients by the same amounts 
in regions in which the heat-transfer coefficients are large. 
Such a region exists at the stagnation point on a cylinder 
with an axis ratio of 1:2, and especially on a cylinder with 
an axis ratio of 1:4. 

The variation in coolant flow required to maintain con- 
stant wall temperature for transpiration-cooled cylinders 
with circular and elliptic cross sections is shown in figure 19. 
The calculations were made for a temperature ratio T;,/T. 
of 1, a value of ¢ of 0.5, and a Prandtl number Pr of 0.7. 
Figure 19 shows that the highest local coolant-flow rates are 
necessary near the stagnation point in order to keep the 
wall temperature down at that place. The magnitude of 
the coolant-flow rate at the stagnation point is proportional 
to the square root of the velocity gradient ди Фи“; this in 
turn is determined mainly by the value of tbe radius of 
curvature at this point. As this radius of curvature de- 
creases, the required coolant flow increases. This is in 
agreement with figure 19, which shows that the maximum 
coolant flow is required at the stagnation point of the elliptic 
cylinder with the 1:4 axis ratio. 
nation point, the flow rates decrease for each cylinder. 
Figure 19 also shows that the use of the thermal rather than 
the convection boundary-layer thickness results in only a 


very minor increase in coolant flow required to maintain the ` 


circular cylinder wall at a constant temperature. 
EXTENSION OF CALCULATION TO HIGH-VELOCITY FLOW 


The heat generated by internal friction was neglected in 
equation (5) according to the assumption of small velocities. 


coefficients meres Дый pan рио сюр elliptic ee 
nis 


Downstream of the stag- |.. 
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The equation 
q—h(T,— Т.) (42) 
gives the heat-transfer coefficient for this case. It was 


already explained that the inclusion of the internal friction 
for a gas with a Prandtl number of 1 results only in the 
change that the temperature T in equation (5) and the 
temperature 7, in equation (42) are now total temperatures, 
as long as the property values may be regarded constant. 
The heat-transfer coefficients determined in this report may 
be used in this case. It is shown in reference 35 by use of 
results obtained m reference 36 that the heat-transfer co- 
efficients determined for low-velocity flow apply to high- 
velocity flow up to a Mach number of about 4 for a gas 
with a Prandtl number different from 1, when the stream 
velocity is constant (flat plate) and the heat flow is not too 
large. The heat-transfer coefficient, however, has now to 
be defined by the equation 


q—h(Taa— Та) 


in which the temperature Tag denotes the value which an 
unheated plate assumes in the high-velocity flow. The 
adiabatic wall temperature may be determined from the 
recovery factor 


(43) 


Е Та- 1, 
v Th T, 


(44) 


which was found to be.equal to Pr for laminar flow and for 
Prandtl] numbers of approximately 1. The difference be- 
tween the total and the static temperatures in the stream 
is connected with the stream velocity by the equation 


(45) 
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FIGURE 19.—Coolant flow required to eer ee cylinder wall temperature. о, 0.5; 


For the flat plate with & constant stream velocity, the adia- 
batic wall temperature is therefore constant. 

Conditions are more involved on a cylinder with a stream 
velocity which varies along its circumference. Even when 
the recovery factor is assumed to be constant, equations 
(44) and (45) give an adiabatic wall temperature which 
varies along the circumference of the cylinder. The fact 
that the low-velocity heat-transfer coefficients also repre- 
sented the high-velocity values on a flat plate, however, 
followed from the fact that the energy equation for constant 
property values is linear in 7, and that a general solution 
of the nonhomogeneous equation describing the heat transfer 
including the internal friction could therefore be obtained 


by superposition of the solution of the homogeneous equation . 


valid for small velocities and a particular solution of the 
nonhomogeneous equation. Such a superposition results in 
a constant wall temperature on the flat plate when the solu- 
tion of the homogeneous equation for constant wall tempera- 
ture and the one describing the adiabatic wall temperature 
is used, since the adiabatic wall temperature is also constant. 
For a cylinder with an arbitrary cross section, however, the 
adiabatic wall temperature which represents a particular 
‘solution of the nonhomogeneous equation varies along the 
circumference. Therefore, a superposition of this particular 
solution with the low-velocity solutions for constant wall 
temperature does not give a constant wall temperature, which 
was specified for the problems investigated in this report. 
Accordingly, the heat transfer has now to be calculated with 


the equation: 
4--(ҚТ, ejf — Tw) 


in which Tar has to be determined for constant wall tem- 
perature conditions; that is, Ter is the temperature which 
& particular spot along the surface, for which the heat- 
transfer coefficlent is to be determined, assumes when the 
heat flow through the wall at this particular spot is zero 
and the wall temperature along the circumference of the 
cylinder is constant. : 

For flow around wedges, this temperature, which may be 
referred to as Ше “effective temperature," can be found 
from the results in reference 23. It is also determined for 
several cases in reference 37. The calculation procedure 
which determines this effective wall temperature from ref- 
erence 23 is described in appendix D. The calculation 
shows that this temperature may be again expressed by a 
recovery factor 


(46) 


Тоу Ly 


ToT р, (47) 


The index « is used to indicate that such a recovery factor 


could be determined experimentally by & model made of a. 


material with a very large heat conductivity so that the 
internal heat conduction would eliminate all temperature 
differences along thesurface. On theother hand, the recovery 
factor describing the adiabatic wall'temperature in equation 
(44) has to be determined experimentally by a model made 
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of & material with an infinitely small heat conductivity so 
as to eliminate internal heat flow. Values for the recovery 
factor r, determining the effective temperature of a wedge 
are presented in figure 20. The recovery factors ro describing 
the adiabatic wall temperature according to equation (44) 
have been calculated for wedges in reference 7. This calcu- 
lation had resulted in values which decreased slightly with 
increasing Euler number m. Repetition of these calculations 
on an electric computing machine, however, according to a 
communication from Arthur N. Tifford of Ohio State 
University, showed that the recovery factors for the adi- 
abatic wall temperature are practically independent of the 
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FIGURE 20.—Recovery factor values determining effective temperature of wedge. 


Euler number and have the same values as the recovery 
factor r, shown in figure 20 for an Euler number m equal 
to zero. 

The consideration up to now dealt with solid surfaces. 
No information was found in the literature on recovery 
factors for transpiration-cooled surfaces. Some recovery 
factors were therefore determined for a transpiration- 
cooled flat plate and a flow with constant property values 
(the same for outside and coolant flow) by an integration 
of the boundary-layer equation (4). The integration was 
carried out in the same way as іп reference 5. The dimension- 
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less stream function f and its second derivative were taken 
from reference 30. The results of this calculation are 
presented in figure 21 and the following table where Т,/Тъ-- 1 
and Pr=0.7: 


Recovery 


The figure shows that the recovery factors decrease consider- 
ably with increasing coolant flow. The calculations were 
extended to positive values of f. which apply to a surface with 
suction. 
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FIGURE 21.—Recovery factors for DU on сово. п plate and constant fluid prop- 
erties. Т/Т.,1; Pr, 0 

It might be worthwhile to mention that the accurate 
determination of the adiabatic or effective wall temperature 
appreciably influences the heat flow as calculated by equation 
(43) only when the difference Т, — T» is of the same order of 
magnitude as or of a smaller order of magnitude than the 
difference Тт,-- T, (see also appendix D). 


RESULTS AND CONCLUSIONS 


Àn approximate method for the calculation of heat transfer 
in the laminar region around cylinders of arbitrary cross 
section was presented. The method, called the equivalent 
wedge-type-flow method, is based on exact solutions. of the 
laminar boundary-layer equations for wedge-type flow and 


245 


takes into account the influence of large temperature differ- 
ences between the flow and the cylinder wall and the influence 
of transpiration cooling. The use of prepared charts reduces 
calculations to a graphical solution of an ordinary first-order 
differential equation. The method can be based either on 
the convection thickness or on the thermal thickness of the 
boundary layer. The results of calculations based on one 
thickness differ slightly from those based on the other 
thickness. ‘There are not enough experimental data avail- 
able to decide which boundary-layer thickness should be 
used. Near the separation point, however, the results ob- 
tained with the thermal boundary-layer thickness seem 
somewhat more plausible. 

. The method was applied to circular and elliptic cylinders, 
дай the following results and conclusions are given: 


1. Results of experiments and exact calculations were 
available only for circular cylinders with solid surfaces. 
Calculations based on the present method and on the thermal 
boundary-layer thickness agreed within 5 percent with the 
exact calculation and within 8 percent with experiment when 
the immediate neighborhood of the separation point was 
excluded. 

2. With.the present method, РСР coefficients may 
be obtained without a knowledge of the flow boundary layer. 
Consequently, such calculations are more rapid than those 
based on the momentum and heat-flow equations. 

3. Heat-transfer coefficients determined from wedge solu- 
tions agreed on the circular cylinder within 15 percent with 
the results of experiments. The kanaka procedure 18 
still more rapid. 

4. For elliptic cylinders, the differences between the results 
of calculations with the various methods decreased as the 
axis ratio increased from 1:2 to 1:4. 

5. The development. of the boundary layer i is determined 
by the velocity distribution around the cylinder. The 
accuracy which has to be expected for the results of calcula- 
tions with the different methods will therefore depend on the 
character of the velocity distribution. 

6. For cylinders with solid walls, the variation of the 
heat-transfer coefficients with ratio of stream to wall tem- 
perature was comparatively small. 

7. For transpiration-cooled surfaces, the effect of tempera- 
ture ratio on heat-transfer coefficients became pronounced, 
especially on cylinders with nearly circular cross sections. 

8. A considerable decrease in heat-transfer coefficients 
accompanied transpiration cooling. 

9. The influence of transpiration cooling on the recovery 
factor was investigated for a flat plate and constant property 
values. It was found that the recovery factor decreased 


considerably with increasing coolant flow. 
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APPENDIX A 
SYMBOLS 


The following symbols are used in this report: 
А dimensionless wall temperature gradient taken from 


| 2 
refs. 22 and 23, тт 


C constant 

» Specific heat at constant pressure 

Е function 

7 dimensionless stream function, (р.у) Nx 

h heat-transfer coefficient 
k 
1; 


с» 


thermal conductivity 
characteristic dimension (major axis of cylinder) 
yt n =M бира, oti*) (see eq. (34)) 


т Euler number, = epee и 
: pss] 


M = 


= От" 


r du; жож 
N =b an =N ( 72% 0*5 vss ) (see eq. (87)) 
Nu Nusselt number, AL/k, 
Pr Prandtl number, c,u/k 
p pressure 
q heat flow 
ў approximated heat flow 
qe heat flow by conduction 
0, heat flow by radiation 
Пе Reynolds number, v, Lp te 
He, Reynolds number, избра ше | 
ry recovery factor defined by (Taa— Т.) Т. Ts) (eq. (44)) 
, recovery factor defined by (Ter — 7,)/(Tr,.—T;) (eq. (47)) 
T temperature in boundary layer 
T, temperature in stream 
и velocity component along surface 
и,  free-stream velocity 
ux dimensionless mass velocity in free stream, р,» 
v velocity component normal to surface 


0% dimensionless velocity normal to surface, 


v 
Ре V Reg 
Potls,o 


7 distance from stagnation point along surface 


z* dimensionless distance from Заря нод point along 


surface, z/.L 
y distance normal to surface 
г dimensionless boundary-layer coordinate taken from 


refs. 22 and 23, 1 | m 7 
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105, [uz 
42—B * 
B pressure gradient parameter, 2m/(m-+1) 
5 boundary-layer thickness 
5% dimensionless boundary-layer thickness, (8/L) Вер 
6, | convection boundary-layer thickness (eq. (15)) 
52 dimensionless convection boundary-layer thickness, 
MURS 
ба | displacement boundary-layer thickness (eq. (13)) 
6, | momentum boundary-layer thickness (eq. (14)) 
5, thermal boundary-layer thickness (eq. (16)) 
ô dimensionless thermal boundary-layer 


(5,/L) УКео 
7 dimensionless boundary-layer coordinate, уч Е 


| Pulls 

6 iba 

ii: Поб 

0 dimensionless temperature-difference ratio, 2 15 
8 к 


2 dimensionless temperature-difference ratio, T dr 
T: 7 
dimensionless stream function taken from ref. 27, 
UTI 
2 w 

absolute viscosity 

kinematic viscosity, u/p 

distance along wedge, taken from refs. 22 and 23 

density 

dimensionless temperature-difference ratio, t 
8 е 


thickness, 


> 


stream function ' 
ubscripts: 

adiabatic 

coolant, when used with 7 

effective 


e LLS 8 DN KR 


оз N’ 
zi 
P. 


except when used with r, refers to & fixed point in the 
stream 
Superscripts: 
m exponent of distance along surface from stagnation 
point for stream velocity, u,— Cz" 
d denotes differentiation with respect to 7 


APPENDIX B 


EVALUATION OF HEAT-FLOW EQUATION 


The energy equation (5) will be integrated along y 
throughout the boundary layer under the conditions of small 
Mach number, constant wall temperature, and constant 
specific heat 


NC 
с 
0 


Sete? Sy uy IN x (Ese 


The first term on the left side can be transformed by partial 


An analogous transformation of the second term ‘and con- 
sideration that the temperature gradient 07/04 is zero out- 
side the boundary layer (for j= о) result in 


az (oul) dy— [^7 9092 ay+ pot)” — 


J| oy m 


The second and fourth terms cancel because of the con- 
tinuity equation (3). Іп the first term, the sequence of 
differentiation and integration can be reversed. Introduc- 
tion of the convection thickness of the boundary layer leads 
finally to the integrated heat-flow equation 


T (0,8) — Putu = ЕЗІ (В1) 


It will now be proved that equation (26), used for the method 
of the equivalent wedge-type flow, is the same as this inte- 
grated heat-flow equation when the convection thickness for 


the boundary layer is used. Equation (Bl) may be trans- 
formed by partial differentiation of the first term into 


ди | ky (90 
sud dz et pide Ge ізді қаты (5% (B2) 


For wedge-type flow, the convection thickness is given by 
the expression 


1-в 
- 2 
ô= А, E ER т 3 
Ne pulls 7 pu (B3) 
Differentiation of this equation gives 
rm 
de 2 "Vp." 9 IN peru, (B4) 


Introducing this expression as well as equations (9) and (12) 
into equation (B2) gives the equation 


ttm К _ Ёш g 
2 (рте Pujs)— p. 0, 


(В5) 


which interconnects the convection thickness with the di- 
mensionless temperature gradient at the wall. The gradient 
of the convection thickness may now be determined from 
the integrated energy equation (B2) when the expressions 
in this-equation are transformed to the new variables . 


dé, | Ho ро Lim 
de my, a] ع2‎ LEM у e 
2 | рег Ps Et Das 28. "o 


Replacing the nondimensional temperature gradient in this 
equation by equation (B5) results in 


dé, 1—m По lcu а Кю 
Pi, 2 iE ридсИз 


which 18 the same ав equation (26). 

It can also be proved by a completely analogous calcula- 
tion that the method of the equivalent wedge-type flow, 
when it is used to calculate the momentum thickness of the 
flow boundary layer, satisfies the integrated momentum 
equation which is obtained from equation (2) by an integra- 


‚ tion over y in a manner similar to the derivation of equation 


(B1) . 
= (р, BOH ри, SE 5, 
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APPENDIX C 
DETERMINATION OF WEDGE SOLUTIONS 


The wedge solutions which were used as a first approxi- | 
mation in figures 11 to 13 can be obtained very easily with 
the use of figure 22 reproduced from reference 9. ‘The hest- 
transfer coefficient has to be determined on a wedge which 
has the same stream velocity and its gradient at the same 
distance from the stagnation point as the real profile. The 
Euler number for this wedge can be found from equation (23). 
In the dimensionless coordinates it is Р 


(C1) 


The parameter fe which determines the coolant flow through 
the porous wall is found from equation (20), which reads, 
when converted to dimensionless quantities, 


___2 „Е 
das m+1 Uo и 


The value (Nu/ A Rey) -/x*/u* can be determined from figure 


(C2) 
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FIGURE 22.—Heat transfer through шшш Бичу layer with and without porous flow 
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22, and Nu/ Re, is finally obtained by multiplication by 
Ми е”. | 

When the temperature ratio ¢ is prescribed, figure 23 
reproduced from reference 9 can be used to obtain the 
parameter f,, for any Euler number m. Equation (C2) then 
determines the value v? and the distribution of the required 
coolant flow along the profile. | 
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FIGURE 23..—Temperature of porous wall (ref. 9). 


APPENDIX D 
DETERMINATION OF EFFECTIVE WALL TEMPERATURE 


It is shown in reference 23 that for high-velocity flow of a 
fluid with constant property values around a wedge with 
constant wall temperature, the temperature field can be 
expressed by the equation 


T= (To— Tr, (1—0 + (Tr ,— J) 9-- T, (D1) 


in which 0 represents the nondimensional temperature field 
for low-velocity flow and 9, the nondimensional temperature 
field for high-velocity flow and a wall temperature equal to the 
total stream temperature. The heat flow from the wall, 
obtained by differentiating equation (D1), is 


qu —Е ез С Tr) (эу Unt) (7). ] 
(D2) 


With the transformations used in reference 28 (see also 
appendix E) 


е‏ ل 
J2—8 V vx‏ 
иа Ох" (D3)‏ 
—.2m |‏ 
m+1‏ 


equation (D2) can be transformed into 


This equation is to be brought into the form 7 
= = (To—T er) (5 (D5) 
A comparison of equations (D4) and (D5) gives 
49 
ТТТ Ta —( r1.) ` (D6) 
dz)» 


321695—55——17 


from which the difference between the total and effective gas 
temperatures can be found. The expression 


d$ 
Prs Тау | 227 C — dz 10 
Ty .—T, =1 on 40 (07) 
dz/, 


defines this temperature difference and the recovery factor 
for the effective wall temperature. The nondimensional 
temperature gradients appearing on the right side of this 
equation are presented in references 23 and 37. In this way, 
the values in figure 20 have been determined. 

To obtain an estimate of the conditions under which the 
difference between the adiabatic wall temperature and the 
effective wall temperature may be neglected, the heat flow 
into the wall will be approximated by the equation 


TT T, ($7). (D8) 


and the error of such an approximation will be determined: 
The ratio of the exact heat-flow equation (D2) to the one 
approximated by equation (D8) is 


T. € — шеке * @—.„) 


Introducing the recovery factor for the adiabatic wall tem- 
perature i | 
Іт 3 — LaF (1—ro) (Tr, s 4) 
gives 
de 1 Era Ta _ 
d. 1-- тт, = * (r& T'o) 
For an Euler number equal to 1, which characterizes flow . 
near a stagnation point and which, according to figure 20, 
shows a large difference between the recovery factors ry and 
ra, the error is smaller than 5 percent when · 
Тоа T 
1. T : Т з 
is larger than 2.5. 
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` APPENDIX E 


~ 


COMPARISON OF VARIABLES 


This appendix gives a comparison of the variables used in 
references 9, 25, and 26 with the ones used in references 22, 
23, and 27. АП of these references deal with wedge-type 
flow. 


The values used in this report are related to the ones in the 
aforementioned references by the following equations: 


зі б. m U£ 
£p ———— — go 


Symbols 
from 
references 
, 22, 23, and 
21 


Relation among 
symbols 
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